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Abstract: We construct an effective Lagrangian for interacting pions with non-relativistic 
00 . energies. The coupling constants can be expressed in terms of the different scattering lengths and 

^ slopes. When used in the calculation of the pionium decay rate, the scattering slope contribution 

^ ' gives a correction of about 8% compared with the lowest order contribution coming from the 

. scattering lengths alone. 
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Q^ ■ Pionium is a hadronic atom of vr^ and vr bound by the Coulomb force. It is highly 

unstable via the strong decay tt^ + 7r~ tt^ + 11^ which probes the low-energy interactions 
Qh! of the pions. As such it can be used to test more accurately the predictions of chiral 

perturbation theory which is an effective theory for QCD at low energies |]l[H2[. It was first 
1^ . constructed by Weinberg who used it at tree-level to calculate the vrvr scattering amplitudes 

""^ ' in agreement with current algebra results Q. Since then the results have been improved 

with one- loop corrections by Gasser and Leutwyler|Q and are now carried to two-loop 
^ ■ order |5|. On the other hand, the experimental values of these scattering amplitudes are 

still very uncertain. For instance, the isospin-zero S'-wave scattering length is known with 
only 20% accuracy]^. 

Recently a lot of interest has been generated by the possibility of a more accurate deter- 
mination of scattering lengths from measurements of the hadronic decay of pionium[0]|^]. 
In order for this to succeed, one must have a complete understanding of the different effects 
acting in the decay process. Since the pions in this hadronic atom are non-relativistic, 
they can be described by an effective theory expanded in terms of operators of increasing 
dimensions involving pion fields and their derivatives. By matching it to relativistic chi- 
ral perturbation theory or experiments, the a priori unknown coupling constants can be 
determined. It will be done in the following to order p^ where p is the momentum of the 
pions. To this order the different scattering amplitudes are characterized by a scattering 
length a and a scattering slope h which then determine the coupling constants. 
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The same, effective Lagrangian also determines tlie dynamics of the bound pions in the 
pionium atom. Since we have a strictly non-relativistic system, we can use the ordinary 
Schrodinger equation to calculate the wavefunctions and there is no need for covariant 
formalisms like the Bethe-Salpeter equation or others. This is in the very spirit of NRQED 
established by Caswell and LepageQ[10| and used with great success for muoniumpl| and 



positronium||l2|. In this non-relativistic framework one can then systematically calculate 
corrections to the different energy levels. In particular, a complex contribution AE signals 
that the corresponding state is unstable with a decay rate given by F = — 2ImA£' and 
thus with lifetime r = l/F. 

The dominant part of the pionium decay comes from the constant part of the ampli- 
tude for TT~^ + '7T~ — > TT*^ + TT*^, 1.6. from the scattering length. In the following we will show 
that the energy dependence of the amplitude, or the scattering slope, gives an additional 
contribution which is around 8% of the leading term. It is an important correction and 
larger than typical electromagnetic corrections which have been considered until now[0]|^]. 

Non-relativistic pions are described by the complex Schrodinger fields tt = (vr^, ttq, 7r_) 
where 7r_|_ annihilates a vr^, vr* creates a vr" and so on. The free fields are described by 
the Lagrangian 

The masses of 7r_|_ and 7r_ are the same and will be denoted by m_|_ while vro has a slightly 
lower mass denoted by mo. We could also include here a relativistic coupling oc vrlV^vrj, 
but we will ignore such small corrections in the following. In the same vain, we will not 
consider electromagnetic effects although they are in general important in the problem 
under consideration. 

For the interacting part Cint we will assume exact isospin invariance and only S'-wave 
interactions. We then find that the lowest order interaction can only involve two possible 
couplings, 

Cinti^r) = Go{7V* ■ 7V){7V* ■ Tt) + Ho{tV* ■ 7V*){7V ■ 7v) (2) 

Thus, we have the full Lagrangian C = Co{'7T^)+Co{TrQ)+Co{TT-)+Cint{'^)- The interaction 
has dimension six and is thus not renormalizable in the ordinary sense. But considered 
as an effective theory, it can be renormalized to every order in the expansion of J^int in 
higher-dimensional operators. It has essentially the same form as a corresponding effective 
theory proposed for non-relativistic nucleons by Weinberg[^] and recently improved by 



Kaplan, Savage and Wise for np scattering[14] and the deuteron|15]. The divergent loop 



integrals can be regularized by a momentum cutoff, but as for most effective theories, it is 
much more efficient to use dimensional regularization with minimal subtraction. We will 
use this method in the following. 

For dimensional reasons we know that the coupling constants Go and Hq must be 
oc 1/m? where the 'heavy mass' m in our case is the pion mass. They can be obtained 
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by matching to relativistic chiral perturbation theory or directly to experiments. Per- 
forming the matching in the first way, we find to lowest order in the expansion of the 
chiral Lagrangian the effective couplings Go = —1/8/^ and Hq = +3/16/^ with the pion 
decay constant /tt = 92.5 MeV. The resulting couplings between pions in different isospin 
channels can now be deduced from Eq.(^) which takes the form 

l^inti'^) = ^^o(7ro7ro7ro7ro) + Bq{t:\tt*_tt+'K^) + -Co(vr;^7rl7ro7ro + 7ro7rQ7r+7r_) 

+ ^Dq{t:\tt\ttj^ttj^ + t:*_ti*_tt^tt^ +2t:\ttIt:j^ttq + 2tt*_ttIt:^'Kq) (3) 

when written out. With the above tree-level values for the two fundamental coupling 
constants, we now have = = Co/3 = —Dq/2 = 1/4/^. 

In order to compare with experiments, we calculate the S-wave scattering amplitude 
r(p) where p is the CM-momentum of the pions. The real part is usually defined by 



Rer(p) = — a + (4) 




in terms of the scattering length a and the slope parameter h which gives the energy 
dependence of the amplitude to lowest order p. With the above values for the coupling 
constants, it is now straightforward to read off Weinberg's scattering lengths in the different 
isospin channels from the Lagrangian Eq.(^). 

Instead of using these results from chiral perturbation theory at tree level, we can 
instead match the coupling constants to the measured cross sections, i.e. to the observed 
S'-wave scattering lengths oo and 02 for isospin 1 = and 1 = 2 respectively and the cor- 
responding scattering slopes. In connection with pionium we will be especially interested 
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Figure 1: Tree-level contribution to the scattering process due to the contact term Cq. 

in the process vr"^ + vr" — > vr*^ + vr*^. To lowest order in perturbation theory the scattering 
amplitude is given by the Feynman diagram in Fig.l which gives 

T(0)(vr+ + vr- ^ vrO + 7r°) = -^(ao - aa) = Co (5) 

We thus have Cq directly expressed in terms of measured scattering lengths. Considering 
related processes, we can similarly obtain the other coupling constants, 

Svr Svr Svr , , 

(00 + 202), 5o = — ^(ao + 02/2), DQ = ^a2 (6) 



3m^ 6m^ 
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when we combine scattering amplitudes with definite isospin. 



So far these relations are only valid at tree level of the effective theory. The scattering 
amplitudes are real and unitarity is thus not satisfied. This can be achieved by going to 
higher orders in perturbation theory. Again considering vr"*" + vr" ^ tt" + tt", we have two 




Figure 2: One-loop contributions to the scattering amplitude. 

one-loop diagrams of the form shown in Fig. 2. One contains vr+vr^ and the other tt^tt^ in 
the intermediate state. Ignoring here the mass differences, they give the correction 

tW (7^+ + ^ + vrO) = - (^BoCo + ^^oCq) /(p) (7) 

where the factor of 1/2 is due to the two identical particles in the tt^tt^ intermediate state. 
The integral over intermediate momenta 

/dPk 1 
(2vr)3 E - kVm. + ie 

where E = p^/mTr is the total CM energy, is seen to be linearly divergent. Using now 



dimensional regularization with minimal subtraction! 14 1, it is simply given by /(p) = 
— i|p|m7r/47r. Using instead a momentum cutoff A, it would contain a term proportional 
with A. This could then be absorbed by renormalization of the coupling constant Cq. 
The net result is either way a purely imaginary result which arises from the unitarity 
requirement, but does not contribute to the scattering length or slope parameter in Eq.(^). 

However, going to two loops as in Fig. 3 we will obtain corrections to the tree-level 
results. The two intermediate bubbles can again contain a 7r'''7r~ or a vr^vr'^ pair. Summing 
up the four contributions from combinations of different bubbles, we then have the next 
order correction 

T(2)(^+ + ^- ^ + vrO) = + ^AoBoCo + ^AlCo + ^o'Co) l\p) 

'ao - a2){al + aoa2 + al)p'^ (9) 
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after regularization. It is seen to give an energy-dependence of the scattering amplitude 
proportional to p'^ and thus contribute to the slope parameter in Eq.(Q). But such an 
energy dependence at two-loop level can also result at tree level from an operator in the 
Lagrangian containing two derivatives. More specifically, for the vr'^'+vr" 7r'^-|-7r'^ channel 
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Figure 3: Two- loop contributions to the scattering amplitude. 



we consider here, we must include higher dimensional operators in the expansion of the 
effective Lagrangian. For S'-wave interactions there is only one such possible operator to 
lowest order in the derivative expansion, 

Cinti-n^-K-^^-K^) = ^C7o(7r;7rl^ovro) + ^Ca^vrlvroV^o + ^'^V* ttovto) + h.c. (10) 

Here the gradient is defined as ^ = l/2(^ — ^). It corresponds to the vertex V2 in Fig. 4 
with the value (p|V2|q) = ^C2(p^ + q'^) in the CM reference frame. Adding this contri- 
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Figure 4: Vertex due to the double-derivative coupling C2. 

bution to the two-loop result Eq.(^) and matching to the definition of the full scattering 
amplitude in Eq.(Q), we have 

6m± m± 



which gives 



2>mt 



(feo - ^2) + (ao - a2)(ao + «oa2 + 02) (H) 



Since this is a coupling constant in the effective Lagrangian, it can also be used for bound 
state problems. There are no 'off-shell' problems in this approach. 

We are now in the position to consider decay of pionium. The ground state with 
wavefunction ^'(r) has the energy E = 2m+(l — a^/8). It will be perturbed by the 
different hadronic contact interactions in Eq.(P). For instance, at tree level we get a real 
energy shift from the elastic coupling Bq. Its magnitude is simply — i?o|^(0)p where the 
wave function at the origin is |^(0)|2 = 73/7^ with 7 = Q!m+/2. It is proportional to the 
scattering length 09-1-02/2. This is the hadronic energy level shift discussed first by Deser, 



Goldberger, Baumann and Thirring|16| 
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At next order in perturbation theory we must evaluate the diagram in Fig. 5 with a 
TT~^TT~ in the intermediate state. It gives also a real, but smaller contribution proportional 
to (oq + 02/2)^. However, the same diagram, but now with a vr'^Tr'' in the intermediate 




TT TT 

Figure 5: One-loop correction to the ground state energy level which gives the decay rate, 
state, is purely imaginary. In the bound state picture it is given by the matrix element 

AE = - f / / -^^*(p) Co — -4^ Co ^(q) (12) 

2 J (2^)3 J (27r)3 J (27r)3 " 2Am - k^/mo + ie " ^ ^ 

as first shown by Labelle|8| using non-relativistic effective field theory. Here Am = — 
rri- gives the energy of the intermediate state and 

*(P) = (^27^ (13) 

is the Fourier transform of the ground-state wavefunction. It gives the probability to find 
the momentum p in this state. Using the regularized value of the integral Eq.(^, we obtain 

AE = -i— Cl \/2Ammo |^(0)P (14) 
Svr 

This imaginary result signals that the ground state is unstable and will decay with a rate 
r = — 2ImAii^ induced by the hadronic coupling Co- With the value given in Eq.(|5|) it 
becomes 



167r 2"^o /2Am o /-,r\ 

*(o)r^\ (^0-02)' (15) 



9m2 y mo 

which is the standard result. 

However, there is some implicit uncertainty here in what value to use for the pion 
mass mjr. It comes from the definition of the scattering lengths for both charged and 
neutral pions. Taking it to be the charged mass m+ for both of them, we can write the 
rate as F = ro(l — 3Am/2m-|-) where 



^ 167r , ^ , m9 / 2Am , ,9 , , 

Since l^m/mj^ = 0.033, the last factor represents a 5% reduction of the main decay rate 
Eq.(|l6|). Such kinematic corrections will be important in a future experimental determi- 
nations of the scattering lengths from the measured pionium decay rate. 
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Evaluating the two-loop correction to the ground-state energy, we find a purely real 
result since each bubble gives an imaginary contribution. Therefore in the order we are 
working at, the two-loop correction to the decay rate is zero. 

The above standard result for the decay rate is due to the constant part of the vr"'' -|- 
7r~ — > vr*^ -|- tt'' amplitude, i.e. the corresponding scattering length. But it has also an 
energy-dependent component parameterized by the scattering slope b in Eq.(|^. We can 
now easily calculate this effect to lowest order in the corresponding derivative coupling C2 
in the Lagrangian Eq.([lO|). It results from evaluating the same diagram in Fig. 5 but with 
one of the Cq vertices replaced with the C2 vertex from Fig. 4. This gives the additional 
contribution 

to the ground-state level shift Eq.([l^). The integrations over momenta k and q are now 
even more divergent than in the first contribution. But again we invoke dimensional 
regularization. We then find that the part coming from the can be neglected since it 
is smaller by a factor a^mo/Am. Thus the integral over q just gives ^'(0). Writing the 
integral over k as 

(fk k2 r d^k 



(27r)3 2Ammo -k"^ + ie J {2tt) 



2Ammo ^ 
2AmmQ — k^ + ie 



(18) 



we see that the last term is zero with dimensional regularization while the first part is just 
the previous integral Eq.(|8|), used in the calculation of the main level shift. In this way 
we obtain the finite result 

AS(2) ^ Q^Q^ /\mmo^2AmmQ |^(0)p (19) 
47r 

Using now Eq.(^) for the coupling constant C2, we thus obtain the corresponding correc- 
tion AF^^) to the decay rate. It can be written as 



Ar(2) Am 
Fo m+ 



«o - 0-2 



+ (Oq + 00^2 + 02) 



(20) 



when we write mo = m-|- to leading order in the mass difference. The experimental 
valuesQ of scattering lengths and slopes are ao = 0.26 ± 0.05, 02 = —0.028 ± 0.012 and 
60 = 0.25 ± 0.03, 62 = —0.082 it 0.008. Using these, we find that the first term is more 
than an order of magnitude larger than the last. Combined, this amount to a 7.6% -|- 
0.4% = 8.0% correction to the main decay rate. On the other hand, using just the tree- 
level valuesQ, we obtain the very similar result 8.6% -|- 0.1% = 8.7%. With values from 
higher order chiral perturbation theory, the result is again not much different. The overall 
hadronic correction is sizable and larger than other known corrections of electromagnetic 
origin § §. 

As a rough check of this rather large correction, we can try to estimate the decay rate 
directly from the matrix element Eq.(^) for vr^ + vr^ — > vr*^ -|- vr*^ with a = (oq — a2)/3 
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and b = (bo — 1)2) /2>. Taking r oc |r|2, we obtain to lowest order in the scattering slope 
the correction factor 1 + 2^ ^^^^ to the standard result Eq.(^). Dividing the average 
momentum ( ) equally between the initial state where it is 7^ and can thus be neglected 
and the final state where it is 2Amm^, we have exactly the dominant term in the more 
accurate result Eq. (pO|) . 
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